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ABSTRACT 

IRAC images of M81 show three distinct morphological constituents: a 
smooth distribution of evolved stars with bulge, disk, and spiral arm components; 
a clumpy distribution of dust emission tracing the spiral arms; and a pointlike 
nuclear source. The bulge stellar colors are consistent with M-type giants, and the 
disk colors are consistent with a slightly younger population. The dust emission 
generally follows the blue and ultraviolet emission, but there are large areas that 
have dust emission without ultraviolet and smaller areas with ultraviolet but 
little dust emission. The former are presumably caused by extinction, and the 
latter may be due to cavities in the gas and dust created by supernova explosions. 
The nucleus appears fainter at 8 fim than expected from ground-based 10 /im 
observations made four years ago. 
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1. Introduction 

Spiral galaxies are valuable laboratories for stellar evolution, stellar population, and 
star formation studies. Detailed observations of the nearest galaxies are especially valuable, 
because they allow us to characterize the integrated properties of the stellar and interstellar 
components and their global relationships, while maintaining sufficient spatial resolution and 
sensitivity to study the properties of individual regions. 

Infrared studies of spiral galaxies yield insight into both the evolved stellar content and 
star formation. Because red giant stars have their peak emission in the near infrared, this 
wavelength range is most directly related to the underlying stellar mass (e.g., Kauffmann 
& Chariot 1998; Bell & de Jong 2001). Emission from the hottest dust grains in star 
forming regions is detectable at wavelengths longer than about 3 /xm, but grain emission 
only becomes a significant fraction of the overall emission at A > 5 fim with the rise of 
the "aromatic features" (Helou et al. 2000; Lu et al. 2003). These features are usually 
attributed to vibrational modes of benzene-ring-based materials often called polycyclic 
aromatic hydrocarbons (PAH — Leger & Puget 1984). This emission is dominated by 
transient heating from single-photon events, so there is no temperature information in the 
emission colors. 

Spitzer Space Telescope observations of spiral galaxies offer a tremendous advance over 
previous capabilities. The vast improvement in sensitivity gives the ability to map entire 
galaxies quickly, and the wavelength coverage allows direct tracing of the stellar and non- 
stellar components (e.g., Pahre et al. 2004, who discuss classification schemes based on these 
components). Spitzer also provides sufficient angular resolution to map individual clouds 
and star forming regions in nearby galaxies; the IRAC resolution of < 2" corresponds to a 
linear scale of <40 pc at a distance of 4 Mpc. 

This paper reports initial observations of M81 with the Infrared Array Camera (IRAC — 
Fazio et al. 2004) on Spitzer. Gordon et al. (2004) report companion observations at longer 
wavelengths, and Helou et al. (2004) report similar observations of NGC 300. M81 is one 
of the nearest large spiral galaxies at a distance of 3.63 Mpc (Karachentsev et al. 2002). 
Its Hubble type is Sab (de Vaucouleurs et al. 1991), and its far infrared luminosity is 
Lfir = 3 X 10^ Lq (Rice et al. 1988). M81 is a member of a group with ~25 companions 
and has recently interacted with M82 and NGC 3077 (Yun, Ho, & Lo 2000) and perhaps 
NGC 2976 (Appleton & van der Hulst 1988). With isophotal major and minor axis diameters 
of 27'xl4', M81 is an excellent target for initial .Spzteer observations. M81 also is the largest 
galaxy in the Spitzer Infrared Nearby Galaxies Survey (SINGS — Kennicutt et al. 2003), and 
these early observations provide a demonstration of the scientific potential of SINGS for its 
full 75-galaxy sample. 
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2. Observations 

IRAC observations of M81 were made on 2003 November 6. Each of 23 positions was 
observed with six 12-s frames: two repeats at each of three dither positions. Positions near 
the center were observed last to minimize the effects of residual images. The visible disk of 
the galaxy was observed by both of IRAC's side-by-side fields of view, but the extended sky 
coverage is toward the northeast for 3.6 and 5.8 /im and toward the southwest for 4.5 and 
8.0 //m. 

The data were processed with Version 8.5.0 of the SSC onhne pipehne, using calibration 
files ('skyflats' and 'skydarks') derived from in-flight data. Further processing was carried 
out in IRAF: 'skydarks' constructed by median-combining data well outside the galaxy were 
subtracted from 3.6 and 5.8/^m frames, a similarly-constructed 'skyflat' was divided into the 

8.0 /xm frames, and background levels of frames in all 4 bands were set to approximately 
zero by subtracting the median sky levels. The SSC "background matching" software was 
used to refine this procedure. Image artifacts^ were corrected (for cosmetic purposes) with 
the IRAF IMEDIT task. Mosaics combining all the data were made with the 2004 January 17 
version of the SSC software mopex. Mosaics for all four bands were made on the same sky 
grid using 0'.'86 pixels. Because cosmic rays have different appearances in the four bands, 
different mopex settings were used to optimize cosmic ray rejection. Figure 1 shows the four 
resulting images.^ 



3. Analysis 

Inspection of Figure 1 shows three distinct morphological components: a smooth 
distribution of starlight prominent in the 3.6 and 4.5 /im images, a clumpy distribution 
of dust emission seen at 5.8 and 8 /im, and a pointlike source at the nucleus seen at 4.5 /im 
and longer wavelengths. These components are discussed in turn below. 

Quantitative analysis was based on fitting elliptical isophotes (Pahre et al. 2004) to 
all four images starting at a radius of 5". The isophotes as a function of semi-major axis 
were then fit to disk -|- bulge models. Similar models were also fit to circular aperture 
magnitudes, and the differences are taken to indicate the uncertainties. Observed quantities 
have no extinction correction; Galactic extinction is neghgible (Schlegel, Finkbeiner, & Davis 



^Artifacts include banding, muxbleed, and column pulldown and are described in the IRAC Data 
Handbook at http://ssc.spitzer.caltech.edu/irac/dh/ . 

^Data files are available via http://cfa-www.harvard/edu/irac/publications . 
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1998), and extinction internal to M81 is patchy and uncertain.^ 



3.1. Stellar Content 

At 3.6 and 4.5 /im, the isophotes are well fit by the disk + bulge models as shown in 
Figure 2. Parameters are given in Table 1. The fitting also determines total magnitudes, 
which are 3.52 and 3.61 (relative to Vega = 0) at isophotal wavelengths of 3.54 and 4.50 //m 
(Fazio et al. 2004). Uncertainties are limited by the Spitzer calihraiion uncertainty of 10%. 
Comparison with the 2MASS Ks total magnitude 3.83 (Jarrett et al. 2003) shows that any 
non-stellar contribution at 3.6 /xm is small. The derived color Ks — [3.6] = 0.31 is the 
same within the uncertainties as K — L = 0.25 ± 0.05 for normal galaxies (Willner et al. 
1984). Likewise, [3.6] — [4.5] = --0.09 is very close to the expected —0.15 for an MO III star 
(M. Cohen, private communication), a spectral type representative of old stellar populations 
(Frogel 1985).^ 

Despite the overall good fit of the [3.6] — [4.5] color to late type stellar colors over most 
of the galaxy, there arc regions that show redder colors. One region is the nucleus, discussed 
below. In addition, the bulge and disk show different overall colors: [3.6] — [4.5] = —0.17 for 
the bulge and —0.04 for the disk. All explanations for this color difference point to a disk 
generally younger than the bulge. A stellar population even slightly younger — dominated 
by K stars instead of M stars — would have [3.6] — [4.5] close to zero. A greater proportion of 
dusty interstellar matter in the disk may also cause reddening or, if the dust is heated, show 
up in emission, but dust is likely to be a small contributor to the color difference.^ There 
is also a small color gradient within the bulge. In contrast to the disk fit, which shows the 
same effective radius at 3.6 and 4.5 /im, the bulge fit shows a larger effective radius at 3.6 
than at 4.5 /xm. This means the inner regions of the bulge are redder than the outer regions 



•^For example, Kaufman, Elmegreen, & Bash (1989b) find Ay > 8 for some of the narrow dust lanes, but 
extinction is much smaller over most of the disk. Ay = 8 corresponds to £'([3.6] — [4.5]) ~ 0.09, which is 
near the limit of detection in the current data. 

^Thc negative, i.e., blue, color for M stars is because of CO absorption in M star atmospheres in the 
4.5 yum band. K stars have weaker CO absorption and are therefore redder than M stars in the [3.6] — [4.5] 
color. 

^If the color difference between disk and bulge is entirely due to emission from transient-heated hot dust 
grains in the disk, and the dust emission follows the oc spectrum found by Helou et al. (2000), the 

dust contribution to the 3.6 /zm emission would be 7% of the disk flux and 4% of the total flux. The Li & 
Drainc (2001) dust model predicts a steeper slope, which would make the 3.6 fim dust contribution 13% of 
the disk emission at that wavelength. 
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and may indicate a younger stellar component near the nucleus. 

The arm/interarm surface brightness ratio is about 0.5 mag at 3.6 /xm, comparable to 
values seen in other galaxies at 2.2 iim (e.g. Rix & Rieke 1993; Rix & Zaritsky 1995; Grauer 
& Rieke 1998). This ratio should be a good measure of the mass density enhancement in 
the arms.^ Nearly the same arm/interarm ratio is found from a blue image (Cheng et al. 
1997) when the exact same areas of the galaxy are compared. This agrees with ratios found 
by Elmegreen, Seiden, & Elmegreen (1989). In fact, the 3.6 /xm contrast is smaller inside 
a 5' radius and larger outside, just like the contrast in blue light. However, the qualitative 
appearance in blue light differs from that at 3.6 /im; the latter shows a smoother distribution, 
being less affected by star clusters or extinction patches. 

3.2. Dust Content 

One of the most striking results from the IRAC images is the complex structure of the 
non-stellar emission, as traced in the 5.8 and 8 fim images. Based on the morphology and 
colors, we assume that the 3.6 fim image is a direct measure of the stellar light. We have 
therefore subtracted the 3.6 /xm image, scaled according to the expected color of an MO star 
(Pahre et al. 2004), from the 5.8 and 8 iim images to produce images of non-stellar emission. 
The 8 iim non-stellar image is shown in Figure 3. Non-stellar flux densities are 2.3 and 5.9 Jy 
at 5.8 and 8.0 //m respectively, 35% and 69% of the emission at these wavelengths. 

Quahtatively the non-stellar emission traces the same spiral structure as other gas 
tracers such as the H I (cf. Rots & Shane 1975; Appleton, Davies, & Stephenson 1981; Adler 
& Westpfahl 1996) and optical dust lanes, as well as young stellar tracers including Ha, 
ultraviolet continuum emission (Marcum et al. 2001), and the radio continuum (Kaufman et 

al. 1989). When the structure of the IR emission is examined in more detail, however, it is 
most strongly correlated with the optical dust lanes (as might be expected). Figure 3 shows 
large regions where there is strong emission at 8 /xm but no detectable ultraviolet emission. 
Detectable dust emission extends to the outer reaches of the star forming disk but not as far 
as the outer H 1 spiral arms (Adler & Westpfahl 1996), showing a decline in the product of 
PAH abundance and intensity of the stellar radiation. 

The non-stellar emission is both more spatially extended than the stellar distribution 



^If short-lived red supergiants contribute significantly to the light in the spiral arms but little to the 
interarm regions, the arm/interarm mass ratio would be smaller than indicated by the relative surface 
brightnesses. However, Rix & Rieke (1993) found for M51 that red supergiants are significant only in a 
single small patch and not in the overall emission. 
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(as shown in Figure 2(d) by the rising ratio of non-stellar to stellar emission beyond 100") 
and much clumpier (as shown in Figure 1). The arm/intcrarm ratio is 3 in the 8 /xm image 
(including starlight), similar to the ratio for NGC 6946 at 7 and 15 /xm (Malhotra et al. 1996). 
The non-stellar image (i.e., with starlight subtracted) is so clumpy that an "arm surface 
brightness" is virtually meaningless, but even modest clumps have surface brightnesses an 
order of magnitude greater than interarm regions,^ and the brightest clumps are more than 
3 times brighter still. 

The interstellar medium emission at the 5.8 and 8 /xm is likely dominated by aromatic 
features usually attributed to PAH, as demonstrated for instance by the ISO-PHT spectra of 
galaxies analyzed by Lu et al. (2003). Rigopoulou et al. (1999) observed the central region of 
M81 (24" X 24" beam) with ISOPHOT-S and detected emission in the 7.7 //m PAH feature 
but found a small feature/continuum ratio of <0.115. The low ratio compared to other 
galaxies may result from the large angular size of M81 relative to the more distant galaxies 
observed and thus the lesser disk contribution in the same fixed beam size. Both starlight 
from the bulge and AGN emission (Section 3.3) in the M81 nucleus will dilute the aromatic 
feature emission, which mostly comes from the disk (Figure 3). Aromatic feature emission 
measured with ISO-CAM in the 6.75 fim filter correlates well with H-alpha emission (Roussel 
et al. 2001) (although the Ha to 15 correlation may be stronger — Sauvage et al. 1996). 
PAH at 6.75 /im was therefore proposed as a reliable star formation estimator for typical 
galaxy disks. However, old, cool stars can also excite aromatic feature emission (Li & Draine 
2002), and such emission is therefore not expected to be a perfect star formation indicator, 
any more than far infrared luminosity or any other indicator. 

Figure 3 shows that the non-stellar dust emission overall tracks the near ultraviolet 
emission (NUV), but there are both regions with dust emission and no NUV and with 
NUV but no dust emission. PAH without NUV can be explained by dust extinction; one 
magnitude of visual extinction, which would heavily attenuate the NUV, corresponds to 
an undetectable 0.01 mag of differential extinction between 3.6 and 4.5 /xm. Seeing NUV 
light but not dust emission is harder to explain, but perhaps cavities cleared by supernova 
explosions are responsible. A signature of such cavities is clumpy Ha (Cheng et al. 1997) 
and 8 /xm emission surrounding an area dark at these wavelengths but bright in NUV. This 
morphology can be explained by a supernova explosion having cleared gas and dust from a 
young star cluster, leaving behind and revealing the hot stars. 



''The interarm regions probably have non-zero surface brightness, but the uncertainties are large and hard 
to quantify because the result depends so critically on what stellar spectrum is subtracted and on the still 
poorly- known IRAC calibration for extended sources. 
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3.3. Nucleus 

The nucleus of M81 has been classified as a LINER (Heckman 1980) or Seyfert 1 (Peimbert 
& Torres-Peimbert 1981) with an estimated black hole mass 7 x 10^ Mq (Devereux et al. 
2003). Like many galactic nuclei with inferred black holes, the M81 nucleus seems relatively 
inactive. Nuclear emission varies by a factor of about two in visible light (Devereux et al. 
2003) and near 10 //m (Grossan, Gorjian, Werner, & Ressler 2001) on time scales estimated 
at about 25 years. The ground-based 10 //m data at ~arcsec resolution suggest a pointlike 
nucleus. 

The new IRAC data suggest that the nuclear source has varied in the previous four 
years. IRAC flux densities in a 3'/5 square aperture (mcZuc/m^' starlight) are 109, 75, 65, and 
56 mJy at 3.6 to 8 jiva respectively with uncertainties of about 10%.^ The 8 /xm flux density 
can be compared with the 1999 10 /xm flux density of 159 mJy (Grossan, Gorjian, Werner, 
& Ressler 2001). The spectral shape of the nuclear emission between 8 and 10 is not 
known, although the ISO spectrum in a 24" beam shows a steep rise (Rigopoulou et al. 1999). 
Nevertheless, it seems hard to account for the difference in measured ffux densities (159 mJy 
at 10 /xm in 1999; 56 mJy at 8 /im in 2003) by spectral shape alone. If the difference is due 
to variability, it will strongly support the idea that we are seeing an AGN component. 

In order to examine the nuclear morphology, the bulge + disk model fits were subtracted 
from the IRAC images at each wavelength. The results are shown in Figure 3. A residual 
pointlike source is quite evident at 4.5 /xm and longer wavelengths and might be present at 
3.6 //m as well, although the brightness of the stellar bulge at 3.6 //m and nearby negative 
values in the difference image make it hard to be sure. Flux densities of the apparent 
pointhke source at the four IRAC wavelengths are 6, 21, 22, and 39 mJy respectively, 
with the first being uncertain by more than a factor of two. Davidge & Courteau (1999) 
found an unresolved source at the M81 nucleus on the basis of J-K color, and Quillen et 
al. (2001) measured a fiux density for the point source of 13 mJy at 1.6 /im. Considering 
the uncertainties and possible variability, the spectral energy distribution from 1.6 to 8 /im 
might be consistent with either a power law [F^ oc i/~°'^) or emission from dust having a 
temperature range of 400-1000 K and emissivity that decreases with wavelength. 

The images also show a filament or bar of material leading from a ringlike structure at 
~1 kpc (1') to an inner disk and then to the AGN. The prominence of the ring at 8 /xm 



^These flux densities are relative to Vega flux densities of 277.5, 179.5, 116.6, and 63.1 Jy and are without 
subtracting any stellar component. They include corrections for IRAC aperture effects and are thus directly 
comparable to ground-based observations made with a large chopper throw. Chopping with a small throw 
(e.g., 9'.'6 as used by Grossan et al.) would give smaller flux densities than with a large throw. 
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suggests that much of its emission is from dust and probably mostly aromatic features. These 
inner features are probably related to the structures seen in HST Ha images (Devereux et al. 
2003). 

4. Conclusions 

1. IRAC images of M81 dehneate three components by both morphology and energy 
distribution: evolved stars organized in bulge and disk components; a clumpy, dusty 
interstellar medium with star forming regions; and a pointlike active nucleus. 

2. The disk is redder than the bulge in the [3.6] — [4.5] color. This could be explained by 
reddening in the disk, by small amounts of dust emission in the disk, or by a younger 
stellar population in the disk. 

3. The arm/interarm surface brightness contrast at 3.6 /xm is 0.5 magnitudes, in agree- 
ment with previous results in blue light. For normal stellar populations, this value 
should represent the stellar mass density contrast. 

4. The nucleus of M81 shows an apparent point source. Its flux density appears to have 
decreased by something like a factor of 3 in the past four years, and its SED can be 
represented by a non-thermal power law or by hot dust emission. 

5. Star- forming regions are identified by their dust emission, probably mostly in the 
aromatic feature bands. Overall the dust emission traces the spiral arms as seen in 
visible and near ultraviolet light, but there arc many regions where dust emission is 
present and NUV is absent and some where the reverse is true. The ubiquity of dust 
emission imphes that absorption affects both the distribution and amount of NUV hght 
seen from star-forming regions. 



This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. Support for this work was provided by NASA through contract 1256790 
issued by JPL/Caltech. Support for the IRAC instrument was provided by NASA through 
contract 960541 issued by JPL. M. A. P. acknowledges NASA/LTSA grant NAG5-10777. 

IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Fig. 1. — IRAC images in the four bands. North is up and east to the left in all panels. 
Grey scale is \/Iog to show faint structure. 
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Fig. 2. — Surface brightness profiles in elliptical isophotes. Solid lines in panel (a) show 
the measured surface brightness profiles for all four IRAC bands. Dashed lines show the 
bulge + disk model for 3.6 fim. Panels (b), (c), and (d) and the solid line in panel (e) show 
the measured surface brightness profiles for three different IRAC colors. The dashed line in 
panel (e) shows the surface brightness profile for the non-stellar emission alone (i.e., after 
starlight has been subtracted) in the [5.8] — [8.0] color. 
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Fig. 3. — Upper panel: 8 ;um non-stellar radiation. Middle panel: composite color image with 
blue from the near ultraviolet (Marcum et al. 2001), green from the 3.6 /im image, and red 
the non-stellar emission from the upper panel. North is 65° clockwise from the top, and the 
images are 20' horizontally. Three regions showing NUV but not 8 fim emission are visible 
in the outer spiral arm above the nucleus. Lower panels: nuclear regions after subtraction 
of the bulge -|- disk model. Images are 5' = 5.3 kpc square, and North is up. Table 1 gives 
model parameters for 3.6 fim; values for other bands are similar. White regions near the 
nucleus represent negative values, i.e., areas where the model surface brightness exceeds the 
observed surface brightness. 
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Table 1. M81 Stellar Properties 



Parameter 


3.6 /xm* 


5^ 


Bulge total magnitude 


4.62±0.05 


8.99 


Bulge half-light semi-major axis qb (arcsec) 


54 


85 


Bulge surface brightness^ (mag arcsec"^) 


14.86±0.03 


20.29 


Bulge ellipticity (1 — ^b/ob) 


0.30 


0.26 


Disk total magnitude 


4.01±0.05 


8.23 


Disk half-light semi-major axis ajj (arcsec) 


242 


288 


Disk surface brightness'^ (mag arcsec^^) 


17.09±0.23 


21.71 


Disk ellipticity (1 — bD/an) 


0.53 


0.53 


Galaxy total magnitude 


3.52±0.02 


7.79 


Bulge to disk ratio 


0.57±0.05 


0.49 



^Indicated uncertainties are those from the fitting procedure; they 
do not include overall calibration uncertainties of approximately 
10%. 

''Parameters from MoUenhoff (2004). Bulge fit at B is based on 
an r~°-^^ profile instead of r~°-^^ as used for 3.6 //m. 



^Mean surface brightness for area within the half-light ellipse. 



